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I.  Overview 

The  objective  of  the  MliRCURY  project  is  the  development  of  a  softvvaie  system  for 
assimilating  meteorological  data  for  a  mesoscale  region  from  a  luimbcr  of  source:;,  ami  fusing 
these  data  as  necessary  to  generate  reliable  values  for  standard  meteorological  parameters  at 
locations  for  which  measurements  are  not  available.  MRRCURY  will  assimilate  data  from  a 
i|uality-assurcd,  leal-time  database  of  measurements  anil  numerical  predictions  maintained  by 
the  Integrated  Meteorological  System  (IMli'l'S)  or  a  similar  system.  MURCURY’s  outputs  will 
include  conventional  three-dimensional  (3D)  grids  and  point  values  of  meteorological  parame¬ 
ters,  as  well  as  i]ualitativc  de.scriptions  of  spatially-located  meteorological  features  such  as 
fiont  boundaries.  'I’hese  outputs  will  serve  as  input  to  a  variety  of  tactical  decision  aids 
(TDAs),  including  TDAs  de.signed  to  aid  a  Staff  Weather  Officer  (SWO).  Use  of  MERCURY 
in  both  prognostic  and  diagnostic  modes  is  anticipated. 

An  initial  MERCURY  protjtype,  MERCURY-1  was  developed  in  1988  in  order  to  refine 
a  requirements  statement  for  MERCURY  that  had  been  developed  from  a  consideration  of  the 
task  environment  (Coombs  et  al.,  1988).  The  design,  implementation,  and  evaluation  of 
MliR(TJRY-l  is  de.scribed  in  U.  S.  Army  Atmospheric  Sciences  Laboiatory  (ASL)  Renort 
('R-88-()().)4-4  (Eield.s,  1988).  This  report  al.so  outlines  a  high-level  design  for  a  i.ccond  proto¬ 
type,  MIiRCURY-2,  to  meet  an  expanded  .set  of  requiiemeiUs  generated  on  the  basis  of  exp<*ri- 
ence  with  MI'-RCURY-l.  'fhe  architecture  of  MERCURY-2  is  de.scribed  in  ASL  Report  CR- 
9()-()().M-I  (Fields,  1990),  and  its  implementation  is  described  in  Fields  et  al.  (1991),  which  is 
included  in  die  Appendix.  The  present  report  describes  the  current  .s.atus  of  MERCURY  (the 
designation  "MERCURY-2"  has  been  di.scontinued),  and  the  method.^  being  used  to  evaluate  its 
performance  The  development  of  cloud  image  analysis  algorithms  for  implementation  in 
MF:R(T)RY  is  de.scribcd  by  Pfeiffer  (1990). 

The  cuiient  MERCURY  implementation  is  being  used  to  piovide  weather  and  terrain  data 
lor  two  iipplication  sy.stems:  the  AIMS  intentions  analy.sis  system  (Coombs  et  al.,  1990),  and 
the  WADIF  cnvironmenial  effects  prediction  .sy.stcm  (Eskridge,  1990).  These  systems  provide 
a  testbed  for  evaluating  the  .software  and  luser  interface  tooks  provided  by  MERCURY. 

II.  Weather  Data  Iiige.stion  and  Archiving 

MERCURY  ingests  weather  data  broadca.st  in  real  time  via  a  .satellite  downlink,  using  the 
Scientilic  Data  Management  (SDM)  software  developed  by  the  UniData  program  of  the 
University  Corporation  for  Atmos|)heiio  Rc,scarch  (Campbell  and  Rew,  1988).  Functions  for 
ingesting  cloud  image  data  in  the  visible  and  near  infrared  channels  obtained  from  the  Geosyn- 
chionous  Operational  EnvironnKntal  Satellites  (GOES)  have  been  added  to  the  .system  in  the 
last  year.  Full-resolution  GOES  images  have  approximately  5  km  re.solution  at  mid-latitudts, 
so  thc.se  ckitii  are  adequate  for  analyzing  cloud  cover  over  mc.soscale  region ,  Analysis  of  these 
images  is  di.scusscd  in  Pfeiffer  (1990). 

A  hardware  and  .software  upgrade  that  will  allow  MI-RCURV  to  ingest  giidded  numerical 
piediciions  generated  by  the  U.S.  National  Weather  .Service,  and  surface  :ind  upiici-air  data 
horn  non-U.S.  stations,  is  in  progress.  These  upgrades  will  allow  MliRCURY  to  be  used  in  a 
piognostic  mode  for  the  Continental  U.S.,  and  to  be  u.sed  diagno.stieally  in  any  areas  of  the 
world  for  which  adeiiutite  weather  and  terrain  data  arc  available. 

An  automated  archiving  system,  which  records  all  data  for  an  indicated  area  (hat  is 
leeeived  during  a  specified  time  interval,  has  akso  been  implemented  in  the  la.st  yeai.  This 
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syslcm  allows  MI2RCURY  to  build  its  own  archival  data  sets.  Weather  data  from  any  source 
can  be  u.scd  as  a  MIIRCURY  archive,  |)rovided  that  it  is  written  in  the  file  format  used  by  the 
MliRCURY  archiver,  which  is  derived  from  the  format  used  in  the  satellite  transmi.ssions 
received  by  MliRCURY.  An  example  record  showing  this  format  is  included  in  the  Appendix. 

III.  (Jcographic  liifonnalion  Sy.sleiii 

MliRCURY  employs  ;i  high-resolution  geographic  infonnation  .system  (CIS)  to  provide 
terrain  elevation  and  land-u.se  dJita  for  the  region  of  interest.  This  system  can  al.so  be  used  to 
maintain  information  about  other  aspects  of  the  terrain,  e.g.  the  pie.sence  of  roads  or  beddings, 
atul  to  locate  objects,  such  as  weather  stations,  with  respect  to  the  terrain. 

A  terrain  elevation  databa.se  for  the  Continental  U.S.  with  30  .second  resolution  in  both 
latitude  and  longitude  (approximately  600  m  x  900  m  resolution  at  mid-latitudes),  which  was 
obtained  from  the  U.S.  Geological  Survey,  has  been  implemented  in  MERCURY  in  the  last 
year.  'I'his  elevation  database  provides  the  underlying  coortlinaie  grid  for  the  GIS  component 
of  MEIU'URY.  lmplementalit)n  of  a  surfaee  cover  dalaba.se,  al.so  obtained  from  the  U.S.  Geo 
logical  .Survey,  is  in  progre.ss.  'I'he.se  daiaba.ses  will  allow  MliRCURY  to  be  tested  at  ii  wide 
range  of  locations,  which  exhibit  a  variety  of  weather  patterns,  in  the  Continental  United 
States.  Digitized  terrain  databases  for  central  Germany  have  also  been  obtained,  and  iiniJle- 
mented  in  the  GIS  (Coombs  et  al.,  1990;  Eskridge,  1990). 

The  MERCURY  u.ser  -  developer  interface  allows  any  data  in  the  GIS  to  be  displayed 
graphically  on  a  con.sole  or  terminal  ntnning  an  X-Windows,  version  11  server.  Functions  for 
iciricving  ;ind  disphiying  data  for  any  location  for  which  data  arc  avadable  have  been  added  to 
the  intciface  this  year.  Functions  for  zooming  and  unzooming  on  a  region,  and  lot  ovei laying 
multiple  data  types  have  also  been  added. 

Graphic  editors  that  allow  the  contents  of  the  terrain  elevation  and  land  u.se  databases  to 
be  altered  by  the  user  have  also  been  implemented  in  the  last  year.  T.hese  editors  allow  data 
for  new  regions  to  be  created;  hence  they  can  be  used  to  develop  terrain  data  sets  for  regions 
for  which  digitized  map  data  are  not  available.  They  al.so  allow  the  data  to  be  edited  as  neces- 
saiy  to  update  features  that  have  changed  since  a  ma|)  was  generated.  The  editors  will  be 
integrated  into  MERCURY  in  the  next  year. 

IV.  Analy.sis  Sysleiii 

The  current  MERCURY  analysis  algorithm  generates  two-dimensional  terrain-following 
grids  of  standard  meteorological  parameters  for  the  region  of  interest,  as  well  as  point  values 
for  standard  parameters,  cloud  cover,  ;ind  present  weather.  The  grid  resolution  is  specified  by 
the  user.  The  analy.sis  algorithm  employed  in  the  current  MERCURY  implementation  is  as 
follows. 

1.  Grid  points  within  at  most  one  grid  .spacing  rom  at  least  one  station  thtit  has  repoited  sui 
face  data  in  the  last  two  hours  are  cla.ssed  as  "noni.solated."  All  other  grid  points  aie 
cla.s.sed  as  "i.solatcd." 

2  Heuristics  re  luscd  to  identify  rcpre.scntativc  reporting  stations,  if  any  arc  available,  for 
the  isolated  stations.  "Rcprescntativcncs.s"  is  defined  in  teims  of  both  distance,  and  simi¬ 
larity  of  loctition  with  respect  to  the  surrounding  terrain.  The  heuristics  being  used  tire 
listed  in  the  Appendix. 
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Data  lor  ihc  isolated  grid  points  for  which  rci)rcscntative  stations  arc  availal)lc  arc 
estimated  from  tlie  data  obtained  at  tlie  representative  stiitions.  Tlie  functions  used  lor 
this  estimation  step  are  listed  in  the  Appendix. 

4.  Conventional  exponentially-weighted  objective  analysis  is  used  to  calculate  values  lor  tlie 
noni.solated  grid  points,  and  for  the  i.solatcd  grid  points  for  which  no  rcprc.scnlaiivc  sia 
lions  could  be  identified.  Both  the  measured  data  and  the  estimated  values  calculated  in 
step  3  are  ti.sed  as  data  values  in  this  calculation. 

5.  Values  for  points  not  on  the  specified  grid  are  estimated  by  exponentially-weighted 
averaging  of  the  values  from  the  four  closest  grid  points. 

().  Pre.sent  weather  values  are  taken  to  be  those  reported  by  the  nearest  reporting  station. 

This  algorithm  can  provide  estimates  for  any  point  in  the  region  of  interest.  The  quality  of  the 
estimate  varies,  depending  on  the  proximity  of  the  point  to  one  or  more  reporting  stations,  the 
local  terrain,  and  the  current  weather  pattern. 

It  is  expected  that  the  perfoimance  of  this  algoiilhm  can  be  im|)rovcd  substantially  by 
ilevelojiing  heuristics  that  u.se  cloud  image  data  to  as.se.ss  repre.scntativenc.ss.  Some  of  the  esti¬ 
mation  functions  can  al.so  be  made  more  realistic.  The.sc  issues  are  discussed  in  Sect.  VI 
below. 

V.  I’erfortnaiice  Evalualioii 

A  self-testing  tilgorithm  has  been  implemented  that  allows  MERCURY  to  be  tested  con¬ 
tinuously,  in  real  time,  in  any  region  for  which  data  are  available.  This  algorithm  is  as  fol¬ 
lows. 

1.  A  region  of  interest  and  grid  resolution  for  testing  are  defined  by  the  user. 

2.  On  every  reporting  cycle  (ciich  lunir  for  the  suiface  sttttion  reports  currently  ingesletl  by 
MliRCURY),  a  grid  is  calculated  for  each  reporting  station  S,  using  all  of  the  data  avail¬ 
able  to  the  system  except  that  reported  by  S.  Values  of  all  parameters  for  the  location  of 
sliilion  S  are  estimated  from  this  grid. 

.f.  liriors  foi  etich  panimeter  are  calculated  as  the  difference  between  the  estimated  value  for 
the  location  of  S  and  the  measured  value  reported  by  S. 

I'his  algorithm  can  be  used  witl  archived  data,  in  which  case  ertor  values  for  an  interval  span¬ 
ning  many  reporting  cycles  can  be  calculated.  The  algorithm  can  al.so  be  easily  modified  to 
assess  the  rate  of  degradation  of  performance  as  d;ita  from  progressively  grctiter  numbers  of 
stations  are  left  out  of  the  grid  and  estimate  calculations. 

'Phis  self-testing  algorithm  has  been  run  for  .sevenil  testing  intervals  in  a  300  x  300  km 
legion  surrounding  the  Los  Angeles  basin  (Fields  et  al.,  1991).  The  Los  Angeles  basin  is  a 
ililTicult  test  case,  since  it  includes  both  complex  mountainous  terrain  and  seacotist.  The  tests 
were  run  both  for  MERCURY,  and  for  an  objective  analysis  procedure  alone.  The  results  of 
the.sc  tests  indicate  that  the  performaiKe  of  MERCURY  is  substantially  better  than  the  perfor¬ 
mance  of  objective  analysis  alone  in  this  region.  The  extent  to  which  MF.RC’URY  ouipeifoims 
objective  analysis  alone  incretises  as  d.ila  liom  additional  stations  are  left  out  ol  the  calcula 
lions,  indicating  that  the  peiformance  degradation  of  MERCURY  is  le.s.s  drastic  than  that  of 
objective  analysis  in  data-sparse  en\  ironments.  Repre.senlalive  data  fiom  these  tests  :ire 
included  in  the  Appendix. 
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Somc  tcsls  have  also  been  pcrfonned  for  the  White  Sands  Missile  Range  (WSMR)  area, 
using  data  collceted  by  Atniospheric  Sciences  Labotalory  personnel.  'I'lie  land  use  d.iia  u.seci  in 
these  tests  are  low  resolution;  howevei,  the  results  of  the  te.sts  show  that  MliRC'URY  inovidcs 
some  improvement  over  objective  analysis  alone. 

Vi.  Kocoinnicndaiions 

The  current  MERCURY  implementation  provides  a  working  basis  for  developing  a 
robust,  general  meteorological  data  interpolation  system.  The  principal  tasks  that  need  to  be 
pursued  tire  as  follows. 

I .  Complete*  implementation  of  high-resolution  digitized  land  use  database  for  the  Continen¬ 
tal  U.S.,  and  for  other  areas  as  data  become  available. 

2  Implement  the  cloud  image  analysis  algorithms  described  by  Pfeiffer  (1990),  tind  integrate 
tlicm  into  the  analysis  .system. 

1.  Improve  the  other  heuristics  and  c.stimation  procedures  ii.sed  in  the  analysis  .system,  using 
testing  to  assess  improvements  in  performance. 

4.  Comptire  the  performance  of  the  system  to  that  of  objective  analysis  alone  loi  multiple 
regions  in  the  Continental  U.S.,  over  the  span  of  a  year  so  that  all  seasons  will  be 
rc|)re.scntcd. 

C’ompletion  of  the.se  tasks  .should  enable  MERCURY  to  perfomi  with  the  accuracy  and  reliabil¬ 
ity  needed  to  provide  adequate  data  to  TDAs  and  other  client  systems. 
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Ahstract:  MERCURY  is  an  integrated  software  system  for  acquiring,  analyzing,  and  extrapolat- 
iiig  meteorological  data  in  real  time  for  regions  a  few  hundred  square  kilometers  in  size. 
MERCURY  employs  a  variety  of  data  structures,  including  real  scalars  and  vectors,  three- 
dimensional  models,  chords,  grids,  and  qualitative  descriptions  to  represent  its  input  data,  and 
both  procedural  and  declarative  methods  to  represent  meteorological  knowledge.  MERCURY 
employs  the  extrapolation  mcthCfd  that  it  is  designed  to  improve  upon,  me.soscalc  objective 
analysis,  when  it  is  likely  to  be  sufficient,  and  supplements  this  method  with  additional  heuris¬ 
tic  and  analytic  methods  when  necessary.  Initial  performance  comparisons  indicate  that 
MERCURY’S  extrapolations  are  better  than  those  obtained  with  objective  analysis  alone  in 
regions  of  complex  terrain  and  surface  cover.  MERCURY  is  implemented  in  C,  with  graphics 
in  X-windows,  and  runs  on  a  Sun  3/60  workstation  under  Unix  4.2  bsd. 
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INTRODUCnON 

Meteorologists  are  often  called  upon  to  solve  a  spatial  extrapolation  problem:  given  the 
weather  conditions  observed  at  one  or  a  small  number  of  locations,  what  are  the  weather  condi¬ 
tions  at  some  other  location  in  the  same  region,  for  which  data  are  not  available?  'fhe  weather 
conditions  of  interest  when  solving  this  problem  are  typically  either  mesoscale  or  smtiller,  i.e. 
have  a  maximum  horizontal  dimension  of  a  hundred  kilometers  or  less  and  a  duration  of  less 
than  a  day  (Orlanski,  1975);  examples  include  snowstorms,  thunderstorms,  fog,  or  strong 
winds.  Such  weather  conditions  may  be  influenced  strongly  by  the  local  terrain;  hence 
knowledge  of  the  terrain,  and  of  its  likely  meteorological  effects,  is  essential  for  successful 
extrapolation.  Even  given  such  knowledge,  however,  the  extrapolation  problem  is  often  very 
difficult,  in  part  because  the  available  data  may  not  be  representative  of  the  location  to  which 
the  extrapolation  is  being  made.  The  difficulty  of  the  problem  is  increased  significantly  when 
it  must  be  solved  for  an  unfamiliar  location,  with  fragmentary  data,  and  under  severe  time  con¬ 
straints,  as  is  typically  the  case  in  emergency  response  situations  such  as  hazardous  materials 
accidents,  and  in  military  situations.  Significant  computer  support  for,  and  if  feasible  complete 
automation  of,  meteorological  data  extrapolation  would  be  very  beneficial  in  such  situations. 

This  paper  describes  MERCURY,  a  proof-of-concept  system  being  developed  to  test 
methods  for  representing  meteorological  and  terrain  data,  and  for  performing  spatial  extrapola¬ 
tion.  The  fundamental  requirements  placed  on  MERCURY  are  that  it  i)  accept  actual  meteoro¬ 
logical  data  as  input,  ii)  perfomi  data  extrapolation  adequately  in  any  geographical  region,  iii) 
be  usable  by  a  nonexpert,  iv)  answer  queries  in  at  most  a  few  minutes,  and  v)  run  on  a  rela¬ 
tively  inexpensive  workstation  (McWilliams  ct  al.,  1989).  The  approach  that  we  have  taken  to 
meeting  these  rci|uiremenis  repre.sents  a  compromise  in  which  accuracy  has  been  traded  for 
Sliced:  MERCURY  is  designed  to  make  informed  guesses,  based  on  heuristic  knowledge,  in 
difficult  cases  instead  of  calculating  an  exact  answer  from  a  detailed  physical  model  of  the 
atmosphere.  MERCURY  does  not,  however,  rely  on  “shallow”  heuristic  knowledge  alone;  it 
also  employs  explicit  geographical  information,  and  relatively  simple  analytical  extrapolation 
methods,  in  situations  in  which  these  can  be  expected  to  yield  useful  results.  This  use  of  mul¬ 
tiple  computational  techniques  requires  a  heterogeneous  architecture  capable  of  supporting  a 
variety  of  data  structures  and  their  associated  operations.  MERCURY  thus  departs 
significantly  from  conventional  meteorological  expert  systems,  which  typically  employ 
declaratively-encoded  heuristics  alone  to  forecast  events  of  a  particular  type,  such  as  thunder- 
sionns  or  fog,  in  a  particular  location  or  type  of  location,  and  which  are  generally  intended  for 
interactive  use  by  an  expert  foreca.ster  (e.g.  Zubrick  and  Riese,  1985;  Swetnam,  et  al.,  1986; 
Elio  and  De  Hann,  1986;  Elio  et  al.,  1987;  Jasperson  et  al.,  1987;  McArthur  et  al.,  1987; 
Robens,  1988;  Dyer,  1989;  Stunder  and  Sletten,  1989;  surveyed  by  Moninger  and  Dyer,  1988; 
Moninger  et  al.,  1989).  MERCURY  is  somewhat  similar  to  the  expert  forecasting  system  of 
Zwack  et  al.  (1989),  which  employs  both  heuristic  and  numerical  techniques  to  predict  low- 
altitude  cloudiness. 

The  next  section  describes  the  task  environment  in  which  MERCURY  operates.  Like 
most  .scientific  data  analysis  task  environments,  the  MERCURY  task  environment  is  open,  in 
the  sense  that  the  system  must  be  able  to  perform  adequately  in  an  arbitrarily  large  set  of  dis 
tinct  situations,  most  of  which  cannot  be  characterized  explicitly  in  a  specification  (Hewitt, 
1985;  Partridge,  1987).  System  performance  in  such  task  environments  can  only  be  judged  as 
adequate  or  inadequate  by  compari.son  with  measurements  of  actual  events;  it  is  not  possible. 
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for  example,  to  construct  artificial  test  cases  for  a  weather  system.  The  third  section  describes 
the  data  and  knowledge  representation  strategies  employed  in  MERCURY,  and  the  tradeoffs 
entailed  by  the  use  of  these  strategies.  MERCURY  employs  latitude/longitude  grids,  real 
scalars  and  vectors,  three-dimensional  geometric  models,  spanning  chords  of  areas,  and  qualita¬ 
tive  descriptions  to  represent  its  input  data,  and  both  procedural  and  declarative  constructs  to 
represent  meteorological  knowledge.  The  fourth  section  describes  the  MERCURY  design,  and 
the  implementation  of  each  of  the  software  modules.  The  fifth  section  discusses  the  perfor¬ 
mance  of  the  current  MERCURY  system.  The  final  section  reviews  the  MERCURY  project 
from  the  perspective  of  the  general  issue  of  designing  intelligent  systems  for  analyzing 
scientific  data. 

THE  TASK  ENVIRONMENT 

Meteorological  phenomena  occur  at  a  variety  of  spatial  and  temporal  scales,  ranging  from 
that  of  long-term  variations  in  global  climate  to  that  of  transient  air  turbulence  tiround  indivi¬ 
dual  surface  features  such  as  buildings  or  trees.  A  rough  classification  of  the  relevant  scales  is 
given  in  Table  1;  Orlanski  (1975)  provides  a  much  more  detailed  classification.  The  weather 
occuiTing  at  particular  locations  in  a  region  may  be  due  to  phenomena  at  any  or  all  of  these 
scales;  one  part  of  a  city,  for  example,  may  experience  an  intense  thunderstorm  during  the  pas¬ 
sage  of  a  low-pressure  system,  while  another  part  of  the  same  city  experiences  only  light  rain. 
Such  local  variations  may  be  due  primarily  to  intrinsic  inhomogeneities  in  the  larger  scale  sys¬ 
tem,  or  to  the  effects  of  the  local  terrain.  Regional  terrain  features  such  as  mountains,  coast¬ 
lines,  deserts,  or  large  urban  areas  can  have  strong  effects  on  local  weather  under  some  larger- 
scale  conditions,  while  having  very  little  influence  under  other  conditions  (e.g.  Retallack, 
1984). 


Scale 

Typical  Horizontal 
Length  Scale 

Duration 

Example 

Microscale 

0.1  km 

Minutes 

Dust  Devil 

Mesoscale 

10  km 

Hours 

Thunderstorm 

Synoptic 

100  km 

Days 

Low-pressure  System 

Table  1:  Scales  of  meteorological  phenomena  relevant  to  the  MERCURY 

project. 

A  meteorologist  responsible  for  a  region  typically  has  access  to  a  fixed  set  of  instruments 
for  measuring  temperature,  wind  velocity,  and  so  forth,  which  are  located  at  particular  posi¬ 
tions  within  the  region,  and  which  obtain  data  for  their  positions  at  particular  times.  Data  for  a 
larger  area  containing  the  region,  e.g.  for  the  country  or  continent,  may  also  be  available.  The 
scale  of  phenomena  that  the  meteorologist  will  be  able  to  observe  is  primarily  determined  by 
the  spatial  and  temporal  resolution  of  these  data.  The  meteorologist’s  task  consists  largely  of 
constructing  a  representation,  either  mentally,  on  paper,  or  using  a  computer  workstation,  of 
the  present  weather  conditions  in  the  region  by  extrapolating  from  the  available  data,  and  then 
generating  predictions  from  this  representation.  The  accuracy  of  the  representation  that  is 
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constructed,  and  of  the  predictions  generated  from  it,  is  strongly  dependent  on  the  representa¬ 
tiveness  as  well  as  the  resolution  of  the  data  that  are  collected;  very  accurate  but  nonrepresen¬ 
tative  data  can  lead  to  a  completely  inaccurate  representation  of  regional  weather.  The  quality 
of  the  representation  also  depends  on  the  accuracy  and  robustness  of  the  available  extrapolation 
methods. 

The  standard  method  for  constructing  a  representation  of  the  weather  pattern  over  a  given 
region  from  a  set  of  measurements  is  a  radially-weighted,  linear  numerical  interpolation  from 
the  available  data  to  a  regular  rectangular  grid  (e.g.  Barnes,  1964).  This  “objective  analysis” 
method  provides  the  baseline  performance  on  which  MERCURY  must  improve  to  be  a  useful 
operational  system.  In  practice,  objective  analysis  has  two  serious  weaknesses,  both  of  which 
are  due  to  nonideal  data  resolution  and  representativeness.  First,  closeness  in  space  is  a  good 
indicator  of  similarity  in  weather  only  in  situations  in  which  the  weather  pattern  varies 
smoothly  as  a  function  of  distance,  and  local  effects  are  relatively  unimportant.  If  the  weather 
pattern  does  not  satisfy  these  conditions,  the  values  of  the  wind  velocity  components,  in  partic- 
uliu",  calculated  for  grid  points  even  relatively  near  one  or  more  measurement  sites  may  be 
grossly  inaccurate.  In  any  case,  the  values  for  grid  points  located  far  from  any  measurement 
sites  are  generally  highly  uncertain.  Second,  most  objective  analysis  procedures  repr^’sent  ter¬ 
rain  at  the  resolution  of  the  interpolation  grid,  usually  simply  as  an  elevation  value  at  each  grid 
point.  This  procedure  effectively  smooths  the  terrain  to  the  re.solution  of  the  grid.  The  inter¬ 
polation  is  thus  blind  to  smaller-scale  terrain  features,  which  may  nonetheless  have  a 
signilicant  effect  on  the  weather  pattern.  These  problems  can,  in  principle,  be  solved  by 
acquiring  data  at  more  sites  and  interpolating  on  a  finer  grid,  as  is  being  attempted  for  example 
in  the  U.S.  Program  for  Regional  Observing  and  Forecasting  Services  (PROFS),  but  this  solu¬ 
tion  is  impractical  for  many  operational  systems. 

It  is  the  responsibility  of  the  meteorologist  to  judge  the  reasonableness  of  an  objective 
analysis  of  a  set  of  measurements,  using  knowledge  of  the  local  terrain  and  large  scale 
meteorological  conditions,  qualitative  observations  of  local  or  regional  conditions,  and  any 
other  available  information.  This  requires  considerable  expertise,  which  is  usually  gained  by 
experience  in  a  particular  locale.  Such  experti.se  is  often  very  region-specific,  and  relatively 
nonportable;  even  experienced  forecasters  may  require  significant  time  to  become  familiar 
enough  with  a  new  location  to  perform  well.  These  problems  are  exacerbated  when  the 
meteorologist  is  a  novice,  is  faced  with  fragmentary  or  otherwise  unreliable  data,  or  is  placed 
under  .severe  time  constraints  in  a  situation,  such  as  an  emergency  response,  in  which  lives  are 
potentially  at  stake.  The  standard  approach  to  this  problem  is  to  provide  the  meteorologist 
with  additional  data,  the  computer  resources  to  compute  analyses  over  finer  grids,  and  interac¬ 
tive  graphic  capabilities  for  merging  or  otherwise  manipulating  displays  of  data  or  analysis 
results.  This  solution  is  embodied  in  a  number  of  milittiry  and  civilian  initiatives,  including 
the  U.S.  Advanced  Weather  Interactive  Processing  System  (AWIPS)  effort  to  develop  a  next- 
generation  workstation  for  operational  forecasters.  We  have  taken  a  different  approach  to  this 
problem  in  the  MERCURY  project,  in  part  because  MERCURY  has  been  specifically  designed 
for  a  military  application  in  which  expertise  is  both  scarce  and  expensive,  and  in  part  because 
all  u.ser-intensive  solutions  face  the  problem  that  additional  data  and  manipulation  methods  can 
accelerate  cognitive  overload,  and  hence  decrease  performance,  especially  in  critical,  stressful 
situations.  Our  goal  has  been,  therefore,  to  develop  a  system  that  automates  as  much  of  the 
expertise-requiring  reasoning  as  possible,  using  the  problem-solving  strategies  employed  by 
meteorologists  as  a  guide  in  the  design  of  efficient  heuristic  extrapolation  methods. 
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The  current  MFiRCURY  system  includes  a  graphic  user  interface  tliat  displays  a  map  ol 
the  region  of  interest,  and  provides  access  to  data  and  analysis  results.  This  interface  is 
intended  both  as  an  aid  for  development  and  testing,  and  as  a  prototype  of  an  interface  suitable 
for  use  by  a  novice  meteorologist  working  in  an  unfamiliar  region,  or  by  a  nonmeteorologisl 
with  some  meteorological  knowledge  who  needs  to  use  meteorological  data,  such  as  an 
emergency-response  planner.  MERCURY  is  not  an  interactive  system  in  the  usual  sense;  it 
queries  its  user  for  values  of  a  set  of  parameters  when  it  begins  an  execution,  but  after  that  it 
ingests  and  analyzes  data,  and  generates  extrapolations  autonomously.  Users  can  query  MER¬ 
CURY  for  extrapolated  data  for  particular  points,  but  cannot  direct  MERCURY’S  problem 
solving  procedures.  This  style  of  interface  avoids  the  delays  and  user  biases  as.sociated  with 
user  interaction  in  conventional  expert  systems,  but  does  not  provide  users  with  any  effective 
control  over  the  system.  Testing  of  MERCURY’S  usability  by  the  intended  users  in  an  opera¬ 
tional  setting  will  be  required  to  evaluate  this  interface. 

DATA  AND  KNOWLEDGE  REPRESENTATION 

Our  choices  of  data  and  knowledge  representations  to  employ  in  MERCURY  have  been 
guided  by  the  input-ouput  requirements  of  the  system,  a  desire  for  the  representations  to  appear 
natural  and  familiar  to  a  meteorologist,  and  the  need  for  generality  and  extensibility.  It  became 
clear  early  in  the  project  that  a  single  knowledge  repre.sentation  strategy  would  not  easily 
satisfy  these  guidelines.  We  decided,  therefore,  to  employ  the  most  natural  representation 
available  for  each  type  of  data  or  knowledge,  and  to  develop  algorithms  as  necessary  to  handle 
multiple  representations. 

A  sample  surface  observation  report  of  the  type  used  by  MERCURY  is  shown  in  Fig.  1. 
The  report  contains  both  quantitative  data,  e.g.  values  of  temperature  and  pressure,  and  quali¬ 
tative  data,  e.g.  estimates  of  cloud  type  and  descriptions  of  present  weather.  MERCURY’S 
task  is  to  produce  reports  of  this  type,  and  similar  reports  of  extrapolated  upper-air  conditions, 
for  locations  at  which  data  are  not  available.  The  need  to  generate  both  quantitative  estimates 
of  real-valued  variables  and  qualitative  estimates  of  discrete,  descriptive  variables  requires  the 
use  of  both  quantitative  and  qualitative  data  representations  and  reasoning  methods.  MER¬ 
CURY  employs  procedurally-encoded  continuous-valued  functions  for  quantitative  reasoning, 
and  declaratively-encoded  heuristic  rules  for  qualitative  reasoning.  Control  is  encoded  declara- 
tively  in  the  analysis  subsy.stem,  which  applies  meteorological  knowledge,  and  procedurally  in 
the  lower-lying  data  acquisition,  management,  and  representation  subsystems. 

Meteorologists  use  many  relatively  simple,  analytical  models,  often  encoded  by  single 
algebraic  formulae,  to  estimate  unknown  values  of  real-valued  variables  from  measured  values. 
'riic.se  models  are  compact  representations  of  meteorological  knowledge  that  can  be  applied  in 
ptirticular  circumstances.  They  are  both  “compiled”  and  “deep”  in  the  senses  defined  by 
Chandrasekaran  and  Mittal  (1983):  knowledge  encoded  by  such  models  can  only  be  accessed 
for  certain  tasks,  but  both  the  models  and  the  tasks  they  support  are  general  enough  that  the 
knowledge  has  wide  applicability.  Models  of  this  type  can  be  applied  directly  to  data  encoded 
as  real  scalars  or  vectors  to  produce  real-valued  results;  the  use  of  such  models  thus  avoids  the 
need  to  represent  continuous  variables  by  arbitrarily  selected  discrete  ranges,  as  must  often  be 
done  in  conventional  expert  systems.  Moreover,  because  such  models  are  used  on  an  everyday 
basis  by  meteorologists,  they  satisfy  the  requirement  of  familiiU'ity  to  the  intended  user  popula¬ 
tion.  A  number  of  these  models  have  been  implemented  in  the  MERCURY  analysis  subsystem 
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(see  below),  either  as  extrapolation  methods  for  particular  variables,  or  to  calculate  values  of 
parameters  for  use  in  qualitative  reasoning. 

Meteorologists  make  extensive  use  of  maps,  both  of  terrain  elevation  and  surface  cover, 
and  of  the  state  of  the  upper  atmosphere.  Weather  prediction  expert  systems  designed  to  work 
in  particular  geographic  regions  typically  represent  the  local  terrain  implicitly,  by  encoding  ter¬ 
rain  effects  directly  in  prediction  heuristics.  This  approach  is  not  practical  in  a  system,  such  as 
MERCURY,  intended  for  use  in  many  distinct  geographic  regions  that  may  have  very  different 
terrain.  A  more  explicit  approach  to  spatial  information  was  taken  by  Elio  and  de  Haan  (1986) 
in  the  design  of  METEOR,  which  employed  labeled,  two-dimensional  (2D)  regions  to  represent 
areas  -  corresponding  to  contours  -  of  significant  cloudiness  or  upper-atmosphere  convective 
activity  (but  not  terrain)  on  a  map.  Heuristic  knowledge  of  storm  potential,  represented  by  a 
set  of  rules,  was  used  to  predict  stomi  severity  from  this  spatial  information.  A  representation 
of  this  type  is  feasible  for  terrain,  but  docs  not  explicitly  include  the  slope  information  needed 
to  calculate  terrain  effects  on,  for  example,  wind  patterns.  In  order  to  provide  a  more  natural, 
higlier-resolution  representation  of  elevation,  slope,  and  aspect,  we  have  tidopted  a  full  three- 
dimensional  (3D)  geometric  representation  for  both  terrain  and  upper-atmo.sphere  features  such 
as  low-pressure  systems.  In  each  case,  the  3D  representation  is  constructed  over  a  gridded 
point  representation  of  the  data  set;  information  about  a  single  point  can,  therefore,  be  retrieved 
without  having  to  traverse  a  3D  model.  Models  of  individual  terrain  or  atmospheric  features 
are  constructed  as  separate,  named  entities  that  can  be  accessed  and  manipulated  as  indepen¬ 
dent  objects.  The  3D  representation  has  the  additional  advantage  of  being  visualizable  in  per¬ 
spective,  allowing  the  user  to  gain  an  intuitive  understanding  of  the  shape  of  both  the  local  ter¬ 
rain  and  the  overlying  upper-air  features. 

The  boundaries  of  areas  overlaid  by  significant  terrain,  and  of  areas  with  particular  surface 
cover,  are  also  represented  by  sets  of  endpoints  of  chords  traversing  the  areas.  Marine,  inland 
lake,  agricultural,  forest,  desert,  suburban,  low-density  urban,  and  high-density  urban  areas  are 
represented  in  this  way.  'I’he  chord  representation  allows  efficient  calculations  of  wliether  a 
point  is  in  a  named  area,  of  the  distance  from  a  point  to  a  named  area,  and  of  the  direction 
from  a  point  to  the  nearest  boundary  segment  of  a  named  area. 

Heuristic  rules  are  employed  in  MERCURY  both  to  select  analytical  models  or  other  rules 
to  execute,  and  to  perfomi  qualitative  extrapolation.  The  selection  rules  encode  infomiation 
about  the  applicability  of  particular  approximations,  models,  or  heuristics,  and  can  be  con¬ 
sidered  metarules  (Clancey  and  Bock,  1988).  They  provide  high-level  declarative  control  in 
the  MERCURY  analysis  subsystem.  The  extrapolation  rules  encode  qualitative  physical 
knowledge.  Only  relatively  general  rules  of  either  type  are  used,  in  order  to  avoid  having  a 
very  large  knowledge  base  of  situation-  or  location-specific  heuristics.  This  strategy  maintains 
generality,  but  at  the  cost  of  having  no  natural  encoding  for  the  local  heuristics  that  experi¬ 
enced  meteorologists  use  in  particular  locations.  We  hope  to  compensate  for  this  lack  by 
implementing  methods  for  using  information  from  additional  data  sources,  such  as  cloud 
images  obtained  from  satellites,  and  improved  analytical  models  in  future  versions  of  MER¬ 
CURY. 
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design  AND  IMPLEMENTATION 

MERCURY  comprises  a  sel  of  subsystems,  which  perfomi  i)  data  acquisition  and 
management,  ii)  feature  identification,  3D  t  lodel  generation,  and  area  rcprcsentatioi',  iii)  data 
analysis,  iv)  user  interface  functions,  and  v)  automated  self-testing.  MERCURY  has  been 
implemented  in  C  for  speed  and  portability;  the  graphic  user  interface  is  based  on  X-windows 
(Scheiller  and  Gettys,  1986).  MERCURY  has  been  tested  on  a  Sun  3/60  workstation  running 
Unix  4.2  b.sd. 

Data  Acquisition  and  Management 

The  data  acquisition  subsystem  of  MERCURY  is  shown  in  Fig.  2.  Surface  and  upper-air 
data  from  the  regularly-reporting  National  Weather  Service  civilian  observation  network  in  the 
continental  U.  S.  are  received  in  real  time,  via  a  communications  satellite  downlink,  from  the 
Domestic  Data  Plus  (DD+)  broadcast  of  the  Zephyr  Weather  Information  Service,  Inc.  of  West- 
borough,  MA.  These  data  are  ingested  in  real  time  and  filtered  from  other,  unused  products  in 
the  DD+  data  stream  by  the  programs  txing  and  txdig,  which  are  components  of  the  Scientific 
Data  Management  (SDM)  system  developed  and  distributed  by  the  University  Corporation  for 
Atmospheric  Research  (UCAR)  Unidata  program  (Campbell  and  Rew,  1988).  txing  and  txdig 
run  continuously  in  background  on  the  workstation  that  serves  as  both  the  data  acquisition  .sy.s- 
tem  and  the  execution  platform  for  MERCURY. 

Following  digestion,  the  data  are  written  into  hourly  surface  and  12-hourly  upper-:iir  data 
flics  by  the  programs  convert  (surface  data)  and  uacvt  (upper-air  data),  which  are  components 
of  the  Weather  Processor  (WXP)  system  developed  by  the  Department  of  Earth  and  Atmos¬ 
pheric  Sciences,  Purdue  University,  and  distributed  by  Unidata,  convert  is  automatically 
invoked  hourly  by  an  execution  shell  that  runs  in  background;  uacvt  is  automatically  invoked 
every  12  hours.  Both  the  surface  and  upper-air  data  are  filtered  to  remove  errors  by  Zephyr 
Weather  Services  prior  to  the  rebroadcast;  no  additional  data  quality  checking  is  done  by 
MERCURY. 

The  converted  surface  and  upper-air  data  files  are  available  for  piocessing  by  the  local 
data  management  software  described  below  as  soon  as  they  are  written,  i.e.  hourly  for  surface 
data,  and  12-hourIy  for  upper-air  data.  The  500  millibar  (mb)  height  and  temperature  data 
present  in  the  upper-air  data  files  are  automatically  gridded  at  230  kib  meter  (km)  horizontal 
re.soluiion  by  the  iipcnlc  program,  which  is  a  component  of  WXP,  whenever  a  new  upper-air 
data  file  is  created.  Additional  data  fields  may  be  gridded;  this  only  requires  altering  the 
automatic  execution  shell  that  controls  upcalc.  These  gridded  data  files  are  also  available  to 
the  local  data  management  software  when  written.  The  gridded  data  are,  in  addition,  automati¬ 
cally  processed  by  the  3D  model  generator  described  below  each  time  a  new  file  of  gridded 
data  is  written. 

'fhe  MERCURY  data  management  subsystem  maintains  and  provides  access  to  both  the 
continuously-updated  meteorological  data  files  and  the  static  terrain  elevation  and  surface  cover 
data  files  used  by  the  MERCURY  analysis  subsystem.  When  MERCURY  is  executed,  the  user 
is  asked  to  provide  latitude  and  longitude  (lat/long)  coordinates  of  a  box  bounding  tl^e  region 
of  interest.  The  data  management  subsystem  reads  a  permanent  data  file  that  specifies  the 
lat/long  locations  of  all  U.S.  weather  stations  from  which  data  are  available,  and  selects  those 
stations  that  are  within  the  region  of  interest.  The  most  recent  data  reported  by  each  of  these 
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stations,  together  with  all  previous  data  from  these  stations  back  to  six  hours  before  the  current 
time,  are  read  from  the  available  meteorological  data  files  into  local  surface  and  upper-air  data 
files.  The  contents  of  these  files  are  continuously  updated  using  a  first-in-first-out  (FIFO)  pro¬ 
tocol  as  new  data  become  available.  The  creation  of  local  data  files  allows  the  maintenance  of 
multiple  reports  for  each  local  station,  and  greatly  increases  the  speed  with  which  local  data 
may  be  retrieved. 

The  only  terrain  database  currently  implemented  as  part  of  MERCURY  contains  1  km  x  1 
km  horizontal  resolution  elevation  data  for  a  1(K)  x  100  km  region  of  the  Los  Angeles,  Califor¬ 
nia  basin  area,  together  with  digitized  surface  cover  data  created  using  a  speciali/cd  editor. 
The  local  data  management  sub.system  cunently  only  provides  access  to  these  data  to  the 
analysis  sub.system;  once  additional  terrain  data  are  available,  data  foi  the  region  of  interest 
will  be  .selected  in  response  to  a  user  specification  of  lat/long  cooidinates  in  the  same  way  that 
the  meteorological  data  are  currently  selected.  Implementation  of  a  digitized  terrain  database 
for  the  continental  U.S.  is  currently  in  progress. 

Feature  Identification,  3D  Model  Generation,  and  Area  Representation 

Point  values  of  meteorological  variables  or  terrain  characteristics,  whether  in  isolation  or 
in  the  context  of  a  regular  grid,  are  of  limited  use  in  much  of  the  qualitative  reasoning  that  is 
used  by  expert  meteorologists  to  understand  the  behavior  of  weather  patterns  in  complex  ter¬ 
rain.  The  locations  and  characteristics  of  identified  features,  such  as  front  boundaries,  moun¬ 
tain  ranges,  or  coastlines,  are  much  more  useljl  for  this  type  of  reasoning.  MERCURY 
employs  a  feature  detection  and  3D  model  generation  subsystem  to  construct  3D  representa¬ 
tions  of  terrain  elevation  and  upper-air  conditions  It  uses  boundaries  as  well  as  the  underlying 
point  data  to  represent  surface-cover  characteristic in  2D. 

The  3D  model  generator  is  based  on  the  Electric  Tinkertoys  (ET)  system  developed  by 
Sodcrlund  and  Pfeiffer  (1989),  which  provides  a  library  of  functions  for  constructing  and  mani¬ 
pulating  data  structures  representing  geometric  entites  (Pfeiffer  and  Soderlund,  1988).  It  is 
used  to  build  3D  geometric  models  of  botli  the  gridded  upper-air  height  fields  produced  by 
iipcalc,  and  the  digitized  terrain  elevation  data.  The  model  generator  includes  a  graphics 
display  facility  ba.sed  on  X-windows,  which  forms  part  of  the  MERCURY  user  interface. 

The  model  generator  is  invoked  automatically  whenever  a  new  gridded  upper-air  data  file 
is  written.  3D  models  i.re  constructed  using  all  grid  points  having  height  values  alvve  or 
below  user  specified  cutoff  values;  this  thresholding  procedure  produces  a  segmentation  of  the 
field  into  individual  high-  or  low-pressct.  features,  respectively,  to  be  represented  as  objects  by 
the  modeler.  The  resulting  models  represent  pressir..  '‘'i.e.'s  .and  troughs,  respectively.  Plots 
showing  typical  height-field  models  both  with  and  w'  ,:  ^ut  segmentation  are  shown  in  Fig  3. 

The  model  generator  is  invoked  to  build  models  of  significant  terrain  whenever  MER¬ 
CURY  is  executed  with  a  new  terrain  data  file.  T  'o  measures,  of  elevation  with  respect  to  sea 
level  and  local  slope,  are  used  to  segment  the  terrain  elevation  data  into  regions  of  significant 
terrain;  cutoff  values  for  these  measures  are  obtained  from  the  user.  The  local  slope  is  calcu¬ 
lated  at  each  grid  point  by  calculating  the  cross  product  of  the  vectors  linking  the  point  to  its 
south  and  east  neighbors;  the  angle  between  this  vector  and  the  horizontal  is  the  local  slope.  A 
grid  point  is  taken  to  be  inside  a  region  of  significant  terrain  if  its  elevation  and  slope  values 
are  both  above  the  user-speCifivd  cutoff  values.  A  plot  showing  significant  terrain  models  gen¬ 
erated  from  the  gridded  elevation  data  for  the  Los  Angeles  basin  is  shown  in  Fig.  4. 
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Bouiidary  representations  are  also  generated  for  each  area  overlaid  by  significant  terrain, 
and  for  each  area  with  identified  surfa.;c  cover,  when  MERCURY  is  executed  with  a  tem.in 
data  file.  A  boundary  is  represented  t}'  the  set  of  pai".  of  endpoints  of  all  of  the  f-torth-South 
and  East-West  chords  spanning  the  area.  This  representation  is  suitable  for  concave  as  well  as 
convex  areas,  and  allows  a  rapid  calculation  of  the  size  of  the  area. 


Data  Analysis 

'I’he  MERCURY  data  analysis  sub.f'.fi.  f. turns  values  of  standard  meteorological 
parameters  -  temperature,  dew  point,  press.'.  r.pscd  and  direction,  and  visibility  -  either 

for  a  point  specified  by  its  lat/long  coord:u..i:s.  ■  >  Cr  all  points  on  a  grid  with  user-specified 
horizontal  resolution  (default  =  25  km  x  2.5  k  o.  A  qualitative  description  of  cloud  type  is 
also  returned  in  response  to  queries  for  point  ..'.ta.  A  hierarchy  of  metarules  is  u.sed  to  deter¬ 
mine  how  these  values  are  to  be  calculated,  u-  '  on  features  of  the  region,  and  of  the  large- 
scale  weather  pattern.  The  basic  decision  im; led  by  these  rules  is  between  situations  in 
which  local  effects  are  likely  to  be  important  n  a«,iermining  regional  weather,  and  situations  in 
which  such  effects  are  likely  to  be  unimportant  If  the  region  contains  either  significant  terra.n 
features,  or  any  marine  area,  the  region  is  classified  as  having  local  effects  regardless  of  .ts 
present  weather.  If  a  region  is  covered  oy  a  low-piessure  system,  or  if  any  station  in  the  region 
reports  severe  weather,  the  region  is  classified  as  having  local  e.ffccts.  Otherwise  the  reg.on  is 
classified  as  free  of  local  effects.  The  mesoscale  objective  analysis  code  OA.SYS  developed  by 
the  U.S.  Army  Atmospheric  Scie.ncos  Laboratory  (Henmi  and  Kirby,  1990)  is  used  to  calculate 
a  grid  of  values  from  the  available  data  in  regions  with  no  local  effects.  OASYS  is  based  in 
part  on  the  KRISSY  objective  analysis  system  developed  by  the  U.S.,  Torest  Service  (bosberg 
and  Sestak,  1986),  which  was  initially  used  with  MERCURY.  A  weighting  is  employed 
for  both  surface  and  upper-air  analyses  in  OASYS.  A  radius  of  influence  can  be  specified,  as 
well  as  a  minimum  number  of  stations  to  affect  each  grid  point.  Upper-air  analyses  arc 
im|)roved  by  .'sing  a  cubic-spline  fit  to  the  measured  upper-air  data  from  each  station;  iliis  pro¬ 
cedure  yields  a  uniform  spreUs.  ‘'points  in  the  vertical  dimension.  Horizontal  interpolation  is 
perfonned  at  whatever  vertical  levels  are  designated  by  the  user.  Only  data  from  stations 
reporting  within  the  last  two  hours  are  included  in  surface  analyses;  stations  reporting  in  the 
last  12  hours  are  used  for  upper-air  analyses.  Queries  for  point  data  at  a  point  P  not  on  the 
grid  are  answered  by  linear,  distance-weighted  interpolation  between  t:ie  data  values  at  all  grid 
points  le.ss  than  or  equal  to  one  grid  spacing  away  from  P.  Because  the  qualitative  cloud  type 
value  is  not  gridded,  this  value  is  always  returned  as  the  same  as  the  value  for  the  nearest 
reporting  station. 


Both  objective  and  heuristic  extrapolation  methods  are  used  in  regions  classified  as  having 
local  effects.  Grid  points  more  than  one  grid  spacing  away  from  any  station  that  has  reported 
data  in  the  last  two  hours  are  first  classified  as  isolated.  MERCURY  estimates  data  for  an  iso¬ 
lated  grid  point  heuristically  if  one  or  more  .stations  are  available  that  are  representative  of  that 


encoded  by  rules,  that  measure  distance  and  direction  to  terrain  features,  coastlines,  and  upper- 
air  weather  features.  Data  from  .st.jiions  repre.sentativc  of  a  point  P  are  corrected,  using  either 
physical  iiKxlels  or  heuristic  factors,  for  differences  between  the  elevation  and  surrounding  sur¬ 
face  cover  at  the  station  and  P.  The  corrcf.cd  data  from  all  stations  rcprc.seutat'.ve  of  P  that  are 
no  irore  than  twice  as  far  from  P  as  the  nearest  such  station  are  averaged,  with  linear  distance 
weighting,  to  yield  the  estimated  data  for  P.  Once  data  for  all  i.soiatcd  points  for  which 
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represeniative  stations  are  available  have  been  calculated,  MERCURY  executes  OASYS  to  cal¬ 
culate  a  surface  objective  analysis,  using  both  the  data  measured  at  stations  reporting  in  the  last 
two  hours,  and  the  estimates  obtained  for  isolated  grid  poin,:.,  as  input  to  calculate  data  for 
botn  the  remaining  isolated  grid  points,  and  for  the  nonisolated  grid  points.  Queries  for  point 
data  are  answered  by  distance-weighted  interpolatu.n  between  grid  points,  as  in  the  previous 
case. 

Because  MERCURY  is  implemented  in  C,  rules  aie  encoded  by  it -then  else  conditionals. 
'Flic  working  memory  of  the  rule  bo.se  is  a  set  of  binary-valued  flag  variables  that  arc  set  and 
pas.scd  between  these  conditionals.  The  result  of  executing  a  rule  is,  in  many  <  ases,  the  execu¬ 
tion  of  one  or  more  numeric..!  functions,  e.g.  for  correcting  temperature  data  ^or  altitude.  The 
distinc'.tion  between  procedural  and  det  'arative  encodings  of  knowledge  and  control  is,  there 
fore,  largely  a  matter  of  intert;retation  by  the  designer  in  MERCURY;  the  distinction  is  not 
reflected  in  anj  explicit  way  m  the  software.  Updating  the  knov  '.edge  base  requires  direct 
modification  of  the  code  in  this  implementation,  but  it  has  the  adv.'’".age  of  a  very  straightfor¬ 
ward  relation  between  the  rules  and  the  numerical  functions  that  the./  contiol. 

OASYS  and  the  smoke  screen  program  obskur,  which  is  a  componer/  of  the  user  inter¬ 
face,  are  the  only  components  of  the  MERCURY  system  not  implemented  in  C;  they  have 
been  left  in  the  original  FORTRAN  and  PASCAL,  respectively.  The  analysis  sub.system  writes 
an  input  file  for  OASYS,  and  calls  OASYS  a.‘‘.  an  external  routine;  it  then  reads  the  output  grid 
file  written  by  OASYS  as  necessity  to  answer  queries  for  point  data.  The  interface  with 
ob.skur  i.s  handled  similarly,  except  that  obskur  is  called  by  the  user  interface,  not  the  analysis 
.sub.system. 

The  analysis  subsystem  is  in  an  early  stage  of  development,  with  only  relatively  simple 
extrapolation  methods  implemented  thus  far.  Our  current  effort  is  focussed  on  testing  and 
incrementally  improving  the  extra|)olation  methotis,  and  on  adding  functions  to  allow  cloud 
images  and  numerical  weather  p  dictions  to  be  used  as  input. 

User  Interface 

User  interaction  with  MERCURY  is  supported  by  a  graphic  user  interface  implemented 
using  X-windows.  This  interface  supports  both  color  and  monochrome  graphics  displays  of  the 
region  of  interest,  provides  access  to  both  meteorological  and  terrain  data  via  interactive  menuv 
and  fill-in  forms,  and  allows  user  control  of  MERCURY  via  a  command  menu.  The  interface 
also  incorporates  a  smoke-screen  prediction  program  obskur,  which  demonstrates  the  capabil¬ 
ity  of  linking  MERCURY  directly  to  a  client  application  program  that  requires  extrapolated 
meteorological  data.  A  typical  interface  screen  is  shown  in  Fig.  5. 

The  interface  allows  the  to  display  and  edit  parameter  values  from  any  measurement 
station,  to  scroll  through  stored  measurements  for  any  .station,  to  create  additional  stations,  and 
to  enter  data  for  the.se  additional  stations.  'Fite  latter  options  are  useful  for  porting  historical 
data  .sets  to  MERCURY  for  testing  purposes.  Data  are  displayed  in  pop-up  fomis,  which  may 
be  .scrolled  and  edited  directly.  Edited  station  data  are  written  to  the  relevant  files  by  the  data 
management  subsystem.  Data  for  stations  created  by  the  user  are  treated  in  the  same  way  as 
real  data  by  the  data  management  component;  however,  such  added  stations  are  marked  to  dis¬ 
tinguish  them  from  real,  and  hence  continuously  updated,  stations. 

The  interface  also  supports  specialized  graphic  editors  for  creating  and  altering  terrain 
elevation  and  land-use  data.  These  editors  require  a  color  or  grey-.scale  display  capability. 
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Both  allovy  editing  of  maps  by  a  painting  process,  in  which  the  mouse  is  used  both  to  select 
elevation  or  land-use  values  from  a  menu,  and  to  assign  these  values  to  either  individual  pixels 
or  to  areas.  The  assigned  values  are  mapped  from  pixels  to  lat/long  coordinates  by  the  data 
management  subsystem. 

The  smoke-screen  prediction  program  obskur  developed  at  the  U.S.  Army  Atmospheric 
Sciences  Laboratory  has  been  partially  ported  to  the  interface  to  demonstrate  the  feasibility  of 
two-way  communication  between  MERCURY  and  client  applications.  When  obskur  is 
invoked  from  a  command  button  on  the  MERCURY  interface,  it  requests  the  user  to  select  the 
lat/long  location  for  which  smoke-screen  prediction  is  to  be  performed,  obskur  sends  the 
lat/long  coordinates  of  the  selected  loc.  tion,  together  with  the  current  time,  to  MERCURY. 
The  MERCURY  analysis  subsystem  calcu...tes  the  values  of  all  meteorological  parameters  for 
the  location  and  time  of  interest,  and  returns  these  to  obskur.  obskur  then  calculates  a 
downwind  smoke  plume  from  these  data,  and  displays  the  result  either  as  a  triangle  represent¬ 
ing  the  predicted  smoke  plume  on  the  map,  or  as  a  text  output  to  the  screen. 

Self  Testing 

MERCURY  includes  a  real-time  self-testing  subsystem,  which  may  be  executed  continu¬ 
ously  as  data  arrive.  This  subsy.‘-<em  tests  MERCURY’S  performance  in  extrapolating  values 
of  surface  meteorological  variable^  to  Ir^ations  for  which  measurements  are  available.  When 
invoked,  the  self-test  system  cycle;'  once  per  hour  through  the  set  of  stations  reporting  data.  At 
each  station  S,  it  executes  the  anal>r.is  subsystem  to  calculate  a  set  of  data  values  for  the  loca¬ 
tion  of  S,  using  all  available  data  except  the  data  from  S.  The  differences  between  the  meas¬ 
ured  and  calculated  data  for  S  are  then  written  to  a  test  log.  Because  self  tests  can  be  run  con¬ 
tinuously,  MERCURY’S  performance  at  night  and  at  other  inconvenient  times  can  be  logged 
without  user  intervention.  However,  the  present  system  ;.ocs  not  provide  adequate  testing  of 
MERCURY’S  performance  in  areas  having  no  stations.  This  can  currently  be  done  only  using 
historical  data  sets  obtained  by  placing  temporary  measurement  stations  in  such  locations. 

PERFORMANCE 

MERCURY’S  performance  in  the  Los  Angeles  basin  region  is  currently  being  evaluated 
by  comparing  the  magnitudes  of  MERCURY’S  extrapolation  errors  with  those  of  OASYS  run 
independently.  Example  results  for  two  stations  in  the  Los  Angeles  region,  tested  at  16:00 
G'l’M,  8  March  1990,  are  shown  in  Table  2.  While  both  systems  suffer  prediction  errors, 
MERCURY  consistently  outperforms  OASYS  in  this  region. 

While  the  performance  results  that  have  been  obtained  to  date  are  from  only  a  single 
region,  and  are  not  representative  of  all  weather  conditions  within  tiiat  region,  they  suggest  that 
MERCURY  is  already  capable  of  some  improvements  over  objective  analysis  alone.  Extended 
testing  in  different  regions  and  sea.sons  will  be  needed  to  determine  whether  this  improvement 
is  statistically  significant,  and  to  rationally  refine  the  nicuiwds  used  by  MERCURY. 
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Station 

Variable 

MERCURY 
.  Error 

OASYS 

Error 

LAX 

Temperature 

6.7% 

63% 

LAX 

Dew  Point 

18% 

67% 

LAX 

Wind  Speed 

0.0% 

>100% 

LAX 

Wind  Dir. 

0.0% 

>100% 

FUL 

Temperature 

9.8% 

28% 

FUL 

Dew  Point 

7.4% 

44% 

FUL 

Wind  Speed 

36% 

49% 

FUL 

Wind  Dir. 

43% 

47% 

Table  2:  Comparison  of  MERCURY  and  OASYS  errors  for  extrapolations  of 
meteorological  data  at  stations  LAX  (Los  Angeles  Airport)  and  FUL  (Fullerton). 

DISCUSSION 

MERCURY  is  an  artificial  intelligence  system  designed  to  solve  a  complex,  but  relatively 
well-circumscribed  scientific  problem.  Such  problems  share  a  number  of  characteristics: 

•  They  are  formally  unspecifiable.  The  problems  for  which  an  AI  system  is  needed  are  pre- 
ci.sely  tlio.se  for  whicli  the  correct  answer  is  unknown,  and  for  which  it  will  not  be  known 
until  the  sy.stcm  has  done  its  work.  The  only  problems  for  which  an  input-output 
specification  can  be  given  in  advance  are  those  for  which  the  system  is  unnecessary. 

•  I'he  system  needs  to  be  able  to  work  with  uninterpreted  physical  data.  The  complexity, 
tedium,  and  time  requirements  of  data  interpretation  are  often  the  major  motivations  for 
automation;  requiring  the  user  to  interpret  the  data  in  order  to  render  it  useful  to  the  sys¬ 
tem  is  self-defeating. 

•  The  data  and  the  required  results  are  very  often  real-valued.  Reducing  the  data  to  a  low- 
resolution  discrete  representation  is  usually  unacceptable. 

•  The  system’s  performance  must  degrade  gracefully.  An  inaccurate  answer,  preferably 
with  some  estimate  of  the  error,  can  nonetheless  be  useful;  a  failure  to  process  certain 
data  altogether  because  they  violate  some  expectation  is  unacceptable. 

•  A  well-defined  strategy  for  evaluating  and  incrementally  improving  the  system’s  perfor¬ 
mance  is  needed.  Such  a  strategy  will  often  require  use  of  experimental  data  as  it  arrives 
for  evaluation. 

The  idealized  knowledge  systems  methodology,  with  its  requirement  that  all  domain 
knowledge  be  encoded  in  an  explicit,  declarative  form  (e.g.  Hayes- Roth  et  al.,  1983)  is  inade¬ 
quate  in  task  environments  with  these  characteristics  (e.g.  Partridge,  1987;  Fields  and  Dietrich, 
1990),  and  expert  .systems  designers  are  routinely  warned  away  from  such  task  environments 
(e.g.  Bobrow  et  al.,  1988).  The  demand  for,  and  the  potential  scientific  and  economic  conse¬ 
quences  of,  such  applications  are,  however,  very  large.  Methods  for  developing  systems  that 
will  perform  adequately  in  task  environments  with  these  characteristics  are  needed. 


-20- 


The  approach  we  have  taken  in  MERCURY  combines  elements  of  standard  scientific 
computing  with  elements  of  knowledge-based  AI.  Viewed  at  the  knowledge  level,  MERCURY 
is  an  automated  tool-user;  its  task  is  to  select  and  apply  an  analysis  method  that  is  appropriate 
to  the  region  and  data  at  hand.  MERCURY’S  strategies  for  applying  the  available  tools  are 
explicitly  patterned  on  those  of  an  expert  meteorologist.  This  approach  lia.*?  the  advantage  that 
.scientific  computing  tools  actually  used  by  meteorologists,  such  as  OASYS,  can  be  incor¬ 
porated  directly  into  the  system;  the  knowledge  that  these  tools  encode  docs  not  have  to  be  re¬ 
represented  in  a  different,  and  quite  possibly  less  appropriate,  notation  to  be  used.  Rules  are 
used  to  encode  strategic  knowledge  of  how  the  available  tools  are  to  be  used,  for  which  the 
knowledge-engineering  methodology  and  the  declarative  representation  are  very  appropriate. 

'file  use  of  geometric  data  structures,  such  as  grids,  3D  models,  and  area  rcpre.scniations, 
is  also  modeled  directly  on  the  use  of  such  representations  by  meteorologists.  These  data 
simcturcs  support  spatial  reasoning  in  a  very  natural  way.  While  the  spatial  reasoning  methods 
implemented  in  MERCURY  thus  far  are  relatively  straightforward,  they  are  similar  in  spirit  to 
the  methods  for  automated  visualization-driven  problem  solving  advocated  by  Abelson  et  al. 
(1989).  As  additional  three  dimensional  information,  such  as  3D  models  of  cloud  masses,  is 
introduced  into  MERCURY,  the  sophistication  of  these  methods  will  increase. 

We  believe  that  heterogeneous  designs  such  as  that  employed  in  MERCURY  will  prove 
useful  for  building  AI  systems  for  other  scientific  data  analysis  tasks.  In  particular,  we  expect 
both  the  explicit  3D  spatial  reasoning  methods  and  the  strategy  of  building  automated  users  of 
available  scientific  computing  tools  demonstrated  in  MERCURY  to  be  generally  useful  in 
building  expert  geographic  information  systems. 
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STATION  ID 
TXPE 

LATITUDE  (NORTH) 
LONGITUDE  (WEST) 
ELEVATION  (METERS) 

TIME  (YYMMDDHH) 
TEMPERATURE  (F) 

DEW  POINT  (F) 

WIND  DIRECTION  (DEG) 

WIND  SPEED  (KNOTS) 
PRESSURE  (MB) 

VISIBILITY  (MILES) 

CLOUD  COVER  (C,  S,  B,  0.) 
CLOUD  BASE  (FT) 

PRESENT  WEATHER 


LAX 

SURFACE 

33.93 

118.40 

30.0 

89082916 

64.0 

58.0 

9.0 

4.0 

113.20 

2.5 

B 

1300.0 

F 


Fig.  1 


Typical  surface  meteorological  station  report 
of  the  type  used  as  input  by  MERCURY. 
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local  point 
dualKise 


3D  500  mb 
aitface  models 


Fig.  2:  Block  diagram  of  the  data  acquisition  system  underlying 

MERCURY.  All  programs  are  executed  automatically  whenever 
appropriate  input  data  are  available. 


L 


j 


(Upper)  Plot  showing  an  unsegmented  model  of  a  500  mb  height  field 
at  230  X  230  km  horizontal  resolution.  (Lower)  Plot  showing  3D 
models  of  low  (left)  and  high  (right)  pressure  features  generated 
from  this  height  field,  using  segmentation  cutoffs  of  (upper  level) 
5800  m  and  (lower  level)  5400  m. 
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Fig.  4:  Plot  showing  3D  significant  terrain  models  for  a  100  x  100  km 

region  of  the  Los  Angeles  basin.  The  elevation  cutoff  value  is 
10  m;  the  slope  cutoff  value  is  5°.  The  San  Gabriel  mountains 
extend  across  the  upper  part  of  the  figure. 
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APPENDIX  B.  ARCHIVE  FORMAT 

WAS  will  write  the  weather  data  to  the  given  output  directory,  but  it  wil 
not  be  in  the  same  format  as  the  converted  input  files.  The  WAS  output 
format  for  the  surface  weather  data  is  the  following: 

yymm.srf  -  file  name 

yy  Last  two  digits  of  the  year, 

mm  Month. 

.srf  Three  letter  extension  which  represents  surface  data. 

The  data  is  group  into  three  lines.  Each  format  is  given  below, 
bine  1 : 

sss  aa.aa  bbb.bb  yymmddhh 

sss  Station  identifier. 


aa.aa 

Latitude 

bbb . bb 

Longitude 

yy 

Last  two  digits 

of  the  year. 

mm 

Month. 

dd 

Day  of  month. 

hh 

Hour  of  weather 

observation  in  GMT. 

Example: 

SDB  34.75  118.73  89040414 


Line  2: 

t.t  d.d  w.w  s.s  p.p  v.v 

t.t  Temperature  in  F.  (-99.0  for  missing) 

d.d  Dewpoint  in  F.  (-99.0  for  missing) 

w.w  Wind  direction  in  degrees.  (-99.0  for  missing) 

s.s  Wind  speed  in  )cnots.  (-99.0  for  missing) 
p.p  Pressure  in  millibars.  (-99.0  for  missing) 
v.v  Visibility  in  miles.  (-99.0  for  missing) 

Example : 

!)7.0  36.0  34.0  10.0  -99.0  -99.0 

Note:  The  numbers  do  not  take  up  a  specific  amount  of  space.  There  will 
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always  be  a  space  between  each  number. 


l.inc  3: 


c  Cloud  cover  (1  letter)  ('  '  for  missing) . 

b.b  Cloud  base  (ceiling)  (-99.0  for  missing) . 

p  Present  weather  (may  be  more  than  one  character)  ('  '  for 

missing) . 

Example : 


Note:  Cloud  base  does  not  take  up  a  specific  amount  of  space.  There  will 
always  be  a  space  on  both  sides  of  the  cloud  base. 

Example  WAS  file  follows: 


LGB  33.81  118.15  90110100 

69.00  53.00  270.00  13.00  1014.80  20.00 

MWS  34.10  118.00  90110100 

50.00  45.00  190.00  7.00  -99.00  -99.00 

LAX  33.93  118.40  90110100 

66.00  56.00  250.00  14.00  1015.40  15.00 

EDW  34.90  117.87  90110100 

69.00  31.00-290.00  10.00  1012.40  25.00 

RIV  33.90  117.25  90110100 

71.00  54.00  300.00  9.00  1014.50  10.00 

EDW  34.90  117.87  90110100 

69.00  31.00  290.00  10.00  1012.40  25.00 

BUR  34.20  118.36  90110100 

68.00  52.00  160.00  10.00  -99.00  -99.00 

CRQ  33.13  117.28  90110100 

-99.00  -99.00  270.00  6.00  -99.00  -99.00 

EMT  34.10  118.03  90110100 

-99.00  -99.00  210.00  7.00  -99.00  -99.00 


ONT  34.05  117.62  90110100 


-SI¬ 


TS. 00  50.00  260.00  14.00  -99.00  -99.00 

PMD  34.63  118.08  90110100 

-99.00  -99.00  240.00  16.00  -99.00  -99.00 

POC  34.10  117.78  90110100 

-99.00  -99.00  220.00  6.00  -99.00  -99.00 

RAI,  33.95  117.45  90110100 

-99.00  -99.00  310.00  15.00  -99.00  -99.00 

TOA  33.80  118.33  90110100 

66.00  56.00  280.00  10.00  -99.00  -99.00 

VNY  34.22  118.48  30110100 

69.00  48.00  140.00  6.00  -99.00  -99.00 

WJF  34.73  118.22  90110100 

67.00  42.00  240.00  16.00  1013.10  40.00 

NZJ  33.67  117.73  90110100 

69.00  53.00  230.00  2.00  1014.10  7.00 

SBU  34.10  117.23  90110100 

71.00  48.00  260.00  8.00  1013.30  20.00 

VCV  34.58  117.38  90110100 

70.00  33.00  180.00  12.00  1011.80  20.00 
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AFFEND'LX  C.  REFRESENTATIVE  STATION  SELECTION 

For  each  isolated  point  P,  we  find  the  representative  stations  S  using 
the  following  heuristics. 

IF  . . . 

m  :  Point  P  and  Station  S  are  both  in  the  same  significant 
terrain  object  or  both  P  and  S  are  outside  of  all  signi:. icant 
terrain  objects. 

AND 

#2  :  There  is  no  point  P'  between  P  and  S  with 
elevation (P' )  >  max (  elevation (P) ,  elevation (S)  ) 

AND 

#3  :  Point  P  and  Station  S  are  both  in  the  same  low  pressure 
object  or  both  p  and  S  are  outside  of  all  low  pressure  objects. 

AND 

H  :  Both  P  and  S  are  within  2  grid  spacings  from  an  area  of 
marine  land  use  or  both  P  and  S  are  more  than  2  grid  spacings  from 
an  area  of  marine  land  use. 

THEN  . . . 


Station  S  is  representative  of  point  P. 
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APPENDIX  D.  ESTIMATION  FUNCTIONS 

The  following  is  the  estimation  functions  used  in  the  MERCURY  analysis 
system  for  isolated  point  P  and  station  S. 


1 .  Changes  due  to  elevation 

The  estimated  temperature  at  P  is 

temperature (P)  =  temperature (S)  + 

(elevation (S)  -  elevation (P) )  /  500m  *  1  degree  C 

The  estimated  pressure  at  P  is 
pressure  (P)  -  pressure  (S)  + 

(elevation  (S)  -  elevation  (P) )  /  8m  *  1  millibar 


2.  Changes  due  to  land  use 

if  landuse(P)  is  one  of  {forest}  and 

landuse(S)  is  one  of  {agricultural,  low_urban,  suburban} 
then 

temperature (P)  =  0.95  *  temperature (S) 

if  landuse(P)  is  one  of  {forest}  and 

.landuse(S)  is  one  of  {high_urban,  desert} 
then 

temperature  (P)  =  0.90  *  temperature (S) 

if  landuse(P)  is  one  of  {agricultural,  low_urban,  suburban}  and 
landuse(S)  is  on<'  of  {high_urban,  desert} 
then 

temperature (F)  -  0.95  *  temperature (S) 

if  landuse(P)  is  one  of  {agricultural,  lov/  urban,  suburban}  and 
]anduse(S)  is  one  of  {forest} 
then 

temperature  (P)  =  1.05  '•  temperature  (S) 

if  landuse(P)  is  one  of  {high_urban,  desert}  and 

landuse(S)  is  one  of  {agricultural,  low_urban,  suburban} 
then 

temperature (P)  =  1.05  *  temperature (S) 

if  landuse(P)  is  one  of  {high__urban,  desert}  and 
landuse(S)  is  one  of  {forest} 
then 
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temperature  (P)  =  1.10  *  temperature (S) 


('Iiaiujea  lot  wiiui  .spet^'d  duo  l.o  laud  ust:  when  P  i i.)ulsid>'  ol  any 
aicjuif leant  terrain  objects. 

if  landuse(P)  is  one  of  {forest}  and 

landuse(S)  is  one  of  {agricultural,  desert,  suburban} 
then 

wind_speed (P)  =  0.50  *  wind_speed(S) 

if  landuse(P)  is  one  of  {low_urban,  high_urban} 

landuse (S)  is  one  of  {agricultural,  desert,  suburban} 
then 

v/ind_speed  (P)  =  0.75  *  wind_speed (S) 

if  landuse (P)  is  one  of  {forest}  and 

landuse (S)  is  one  of  {low_urban,  high_urban} 
then 

wi nd_speed (P)  =  0.75  *  wind_speed (S) 

if  landuse (P)  is  one  of  {low_urban,  high_urban} 
landuse (S)  is  one  of  {forest} 
then 

v/ind_speed  (P)  =  1.25  *  wind_speed  (S) 

if  landuse (P)  is  one  of  {agricultural,  desert,  suburban}  and 
landuse (S)  is  one  of  {forest} 
then 

v;ind_speed  (P)  =  1.50  *  wind_speed(S) 

if  landuse  (P)  is  one  of  {agricultural,  desert,  suburban}  and 
landuse  (S)  is  one  of  {low_urban,  high_urban} 
then 

v/ind_speed  (P)  =  1.25  *  v/ind_speed(S) 


•1 ,  Changes  in  v;ind  speed  and  v;ind  direction  when  P  is  inside  a 
significant  terrain  object. 

[A]  if  no  upper  air  data  are  available  for  the  region  and  S  is  lov;er  in 
altitude  than  P,  estimate  wind  direction  at  P  is  eastward  (Northern 
Hemisphere)  or  westward  (Southern  Hemisphere) .  Estimate  v/ind  speed  as 
5  m/s  per  500m  altitude  difference,  added  as  a  vector  to  value  from  S. 

[B}  if  upper  air  data  are  available  for  the  region,  calculate  wiri  speed 
and  direction  at  P  as  functions  of  altitude  from  sounding  data.  If 
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more  than  one  sounding  are  available,  do  for  each  sounding,  then  do 
distance  weighted  average  of  the  results  based  on  distances  to  the 
sounding  stations. 

if  Innduse(P)  is  forest 
then 

wind  speed (P)  =  0.75  *  (estimate  from  [A]  or  [13]) 
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APFENDIX  E.  lA  BASIM  TEST  RESULTS 


Los  Angeles  Basin  Self  Test  Results 


There  are  two  columns  of  numbers  under  each  category.  The  first  is 
the  observed  reading  at  the  NWS  weather  station.  The  second  is  the 
error  between  the  observed  reading  and  the  analysis  systems  output. 

For  example,  in  the  first  row  of  data,  the  temperature  at  LAX  was  75.0 
F.  The  error  between  75.0  F  and  the  output  from  Mercury's  analysis 
systojm  is  -2.3  F.  Therefore,  the  output  from  the  analysis  system  was 
72.7  F.  A  blank  in  the  data  column  indicates  missing  data.  A  bJank 
in  the  results  columns  only  indicates  a  location  for  which  the  tested 
version  of  MERCURY  produced  no  result,  due  to  insufficient  input. 


Station  Date  Temperature  Dewpoint  Hor.  Wind  Ver.  Viind 


LAX 

90071616 

75.0 

-2.3 

62.0 

LGB 

90071616 

72.0 

0.6 

64.0 

NUC 

90071616 

NZJ 

90071616 

72.0 

1.4 

61.0 

RIV 

90071616 

75.0 

-1.0 

63.0 

SBD 

90071616 

73.0 

63.0 

VCV 

90071616 

77.0 

56.0 

BUR 

90071616 

73.0 

-0.5 

64.0 

CNO 

90071616 

CRQ 

90071616 

EMT 

90071616 

FUL 

90071616 

74.0 

-1.3 

63.0 

ilUR 

90071616 

PMD 

90071616 

POC 

90071616 

RAI. 

90071616 

SNA 

90071616 

73.0 

-0.5 

64.0 

TOA 

90071616 

70.0 

3.6 

VNY 

90071616 

72.0 

1.2 

63.0 

EDW 

90071616 

79.0 

59.0 

ONT 

90071616 

77.0 

-3.5 

60.0 

SMO 

90071616 

74.0 

-0.3 

64.0 

WJF 

90071616 

78.0 

1.0 

52.0 

NUC 

90071617 

LGB 

90071617 

73.0 

1.0 

65.0 

LAX 

90071617 

76.0 

-1.0 

63.0 

NZJ 

90071617 

75.0 

0.2 

62.0 

RIV 

90071617 

79.0 

-0.5 

62.0 

SBD 

90071617 

78.0 

64.0 

BUR 

90071617 

76.0 

-0.1 

63.0 

CNO 

90071617 

1.9 

-4.0 

2.3 

6.9 

-5.6 

-1.0 

1.6 

-3.9 

8.9 

-4.7 

2.6 

2.6 

-1.9 

3.1 

2.2 

-0.8 

1.5 

-2.0 

-1.3 

1.4 

-0.0 

-0.0 

-0.0 

-0.0 

• 

o 

1 

0.0 

-1.3 

4.0 

-2.7 

-0.0 

-1.8 

-0.0 

1  .3 

5.2 

-8.4 

3.0 

2.3 

-0.0 

0.4 

-0.0 

2.1 

0.1 

1.2 

-1.6 

6.9 

-2.9 

-0.0 

-4.0 

-0.0 

5.9 

-4.0 

0.3 

6.9 

-4.6 

-0.0 

-1.0 

-0.0 

1  .8 

-0.0 

-0.3 

-0.0 

0.4 

-2.3 

1.2 

-0.1 

6.9 

-2.8 

-8.7 

7.8 

5.0 

-0.8 

0.8 

-0.0 

-1.3 

-0.0 

3.5 

0.9 

-4.9 

4.9 

2.0 

2.8 

-3.8 

3.9 

1.4 

-1.0 

-1.4 

-4.7 

2 .5 

1.7 

2.5 

7.0 

-4.7 

1.5 

1.7 

4.9 

-1.6 

1.4 

-2.8 

7.9 

-2.5 

1.5 

-3.5 

3.0 

6.1 

-0.1 

2.4 

0.0 

1.6 

4.0 

3.1 

1.0 

1.7 

-1.1 

-1.0 

1.5 

-0.0 

-0.0 

-0.2 

0.0 

-5.0 

5.0 

-0.2 

-0.0 

4.0 

-0.0 

3.1 
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CRQ 

90071617 

EMT 

90071617 

FUL 

90071617 

75.0 

-0.6 

63.0 

HHR 

90071617 

ONT 

90071617 

80.0 

-3.4 

60.0 

PMD 

90071617 

POC 

90071617 

RAL 

90071617 

SMO 

90071617 

77.0 

-1.6 

65.0 

SNA 

90071617 

75.0 

-0.3 

65.0 

TOA 

90071617 

71.0 

4.0 

VNY 

90071617 

76.0 

-0.0 

62.0 

EDW 

90071617 

83.0 

59.0 

VCV 

90071617 

81.0 

58.0 

WJF 

90071617 

LGB 

90071618 

75.0 

0.4 

64.0 

LAX 

90071618 

76.0 

-0.2 

63.0 

RIV 

90071618 

84.0 

-0.7 

62.0 

EDW 

90071618 

88.0 

57.0 

BUR 

90071618 

CNO 

90071618 

CRQ 

90071618 

EMT 

90071618 

FUL 

90071618 

76.0 

0.5 

63.0 

HUR 

90071618 

ONT 

90071618 

84.0 

-4.5 

60.0 

PMD 

90071618 

POC 

90071618 

RAL 

90071618 

SMO 

90071618 

77.0 

-1.1 

65.0 

SNA 

90071618 

76.0 

3.0 

64,0 

TOA 

90071618 

73.0 

2.9 

VNY 

90071618 

78.0 

-1.8 

63.0 

NUC 

90071618 

NZJ 

90071618 

81.0 

-4.5 

63.0 

SIM> 

90071613 

83.0 

64.0 

VCV 

90071618 

84.0 

57.0 

MWS 

90071618 

76.0 

1.1 

54.0 

VWF 

90071618 

NUC 

90071619 

NZJ 

90071619 

81.0 

-2.4 

63.0 

VCV 

90071619 

84.0 

57.0 

MWS 

90071619 

76.0 

3.9 

54.0 

LAX 

90071619 

77.0 

0.4 

63.0 

LGB 

90071619 

76.0 

1.4 

64.0 

EDW 

90071619 

93.0 

56.0 

RIV 

90071619 

87.0 

0.5 

62.0 

5.0 

-6.2 

0.0 

7.1 

3.8 

0.7 

1.4 

0.5 

1.0 

-1.2 

2.6 

6.9 

-1.9 

-0.9 

-1.8 

4.9 

1.2 

3.1 

5.4 

-3.8 

4.5 

-4.0 

-1.4 

1.4 

3.8 

-3.8 

8.0 

-5.5 

0.0 

3.2 

-0.0 

2.1 

-0.0 

1.2 

-2.0 

3.8 

-6.6 

10.3 

-5.1 

-2.5 

1.7 

-1.5 

9.8 

-5.3 

-6.1 

5.4 

3.5 

2.8 

1.4 

-9.4 

9.6 

3.4 

2.3 

-0.0 

-1.4 

-0.0 

3.8 

-0.0 

-0.0 

-1.6 

0.0 

-0.5 

5.0 

■3.0 

1.1 

11.8 

-11.8 

2.1 

0.0 

1.0 

4.3 

1.0 

-2.5 

1.2 

2.6 

-3.9 

3.1 

0.7 

-0.0 

4.7 

-0.0 

2.7 

4.5 

-3.1 

5.4 

-1.9 

2.5 

-1.0 

4.3 

-1  .8 

-1.1 

-0.0 

1.2 

-0.0 

3.4 

-0.0 

6.3 

-0.0 

2.9 

1.3 

5.1 

-3.2 

6.1 

-6.7 

5.9 

-3.3 

-1.0 

4.1 

4.3 

-1.8 

-2.5 

5.5 

7.8 

-5.0 

-4.5 

5.2 

-2.4 

4.0 

0.8 

6.9 

-4.3 

-1.0 

0.0 

-0.4 

10.0 

-8.2 

-8.9 

12.0 

1.6 

1  .2 

-0.1 

0.0 

3.3 

10.0 

-6 . 1 

0.6 

-0.7 

1.7 

1.9 

4.0 

4.0 

0.0 

0.8 

-0.6 

8.9 

1.5 

1.0 

2.6 

1.5 

1.5 

-0.7 

0.8 

3.9 

1.9 

-0.6 

5.8 

9.2 

1.5 

2.2 

2.6 

-1.1 

1.0 

7.5 

-7.7 

2.7 

-0.2 

-1.1 

2.7 

-3.5 

7.5 

-4.0 

3.9 

-3.9 

4.6 

-4.6 

2.0 

6.9 

-0.5 

-4.0 

2.5 
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SBD  90071619 

87.0 

65.0 

BUR  90071619 

84.0 

-4.7 

63.0 

CNO  90071619 
CRQ  90071619 
EMT  90071619 
FUL  90071619 

78.0 

-0.1 

64.0 

HHR  90071619 
ONT  90071619 

89.0 

-7.6 

61.0 

PMD  90071619 
POC  90071619 
RAL  90071619 
SMO  90071619 

77.0 

1.3 

65.0 

SNA  90071619 

78.0 

1.9 

65.0 

TOA  90071619 

76.0 

1.2 

VNY  90071619 

81.0 

0.7 

62.0 

WJF  90071619 

90.0 

3.0 

50.0 

MWS  90071620 

76.0 

4.8 

54.0 

LGB  90071620 

78.0 

-0.8 

64.0 

LAX  90071620 

76.0 

2.0 

63.0 

NZJ  90071620 

79.0 

0.1 

63.0 

RIV  90071620 

90.0 

-1.2 

63.0 

SBD  90071620 

88.0 

65.0 

VCV  90071620 

92.0 

54.0 

BUR  90071620 

84.0 

-2.7 

63.0 

CNO  90071620 
CRQ  90071620 
EMT  90071620 
FUL  90071620 

79.0 

-0.7 

63.0 

HHR  90071620 
ONT  90071620 

91.0 

-9.0 

62.0 

PMD  90071620 
POC  90071620 
RAL  90071620 
SMO  90071620 

78.0 

0.2 

65.0 

SNA  90071620 

78.0 

1 .2 

65.0 

TOA  90071620 

75.0 

2.8 

VNY  90071620 

84.0 

-2.3 

62.0 

EDW  90071620 

93.0 

54.0 

NUC  90071620 
WJF  90071620 

93.0 

-3.9 

52.0 

SDB  90071620 

83.0 

56.0 

LGB  90071621 

78.0 

0.8 

64.0 

LAX  90071621 

77.0 

1.4 

64.0 

NZJ  90071621 

80.0 

-0.7 

64.0 

RIV  90071621 

92.0 

0.0 

65.0 

SBD  90071621 

92.0 

66.0 

VCV  90071621 

92.0 

53.0 

4.9 

-0.9 

-0.9 

-3.8 

4.9 

3.2 

-1.2 

6.9 

-0.6 

4.0 

-4.1 

9.8 

-8.9 

1.7 

3.2 

3.8 

-2.2 

1.4 

1  .2 

-1.9 

0.7 

1.8 

3.9 

0.2 

-0.0 

3.6 

-0.0 

3.7 

O 

• 

00 

8.0 

-2.2 

0.0 

2.0 

-0.0 

1.5 

-0.0 

-1.4 

5.0 

-0.1 

0.0 

2.0 

6.9 

-0.7 

-4.0 

4.4 

-2.4 

3.1 

-0.8 

8.5 

-6.3 

-1.8 

2.1 

-1.6 

11.8 

-8.9 

-9.4 

12.0 

3.4 

0.1 

1.0 

-0.0 

-1.4 

-0.0 

3.7 

6.0 

1.0 

-0.4 

-2.8 

3.6 

9.1 

1.5 

2.4 

2.6 

0.6 

-1.4 

4.6 

-3.2 

3.9 

0.0 

0.9 

8.5 

-2.9 

3.1 

-0.6 

1.4 

8.0 

-4.2 

0.0 

7.9 

1.2 

6.9 

-0.3 

-1.2 

-0.4 

2.0 

0.0 

1.7 

-1.0 

-0.9 

-3.0 

5.2 

5.2 

-3.2 

8.7 

0.7 

5.0 

-2.8 

6.1 

-5.1 

-3.5 

8.4 

3.8 

-2.2 

3.2 

-0.6 

-0.8 

0.0 

4.7 

6.0 

-2.4 

9.0 

-4.3 

0.0 

3.7 

-0.2 

13.2 

-6.3 

4.8 

-2.4 

5.2 

-0.4 

-3.0 

2.2 

5.0 

2.0 

0.0 

3.7 

7.7 

0.2 

-6.4 

8.9 

-2.4 

5.0 

-0.2 

8.7 

-6.2 

-1.9 

2.1 

3.5 

11.8 

-9.7 

-7.5 

13.7 

2.7 

0.3 

1.0 

-0.0 

0.0 

-0.0 

5.0 

-0.0 

5.2 

-0.0 

-3.0 

2.8 

5.9 

-2.5 

-1.0 

0.2 

-5.1 

14.1 

-0.2 

0.0 

6.2 

8.0 

-6.0 

0.5 

6.4 

0.7 

7.7 

-5.5 

0.4 

5.6 

-1.4 

-2.1 

8.4 

-0.0 

7.9 

0.6 

-1.4 

0.7 

4.7 

1.7 

-0.7 

1.9 
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BUR 

90071621 

85.0 

-2.7 

64.0 

-0.5 

-1.4 

2.8 

7.9 

-1.8 

CNO 

90071621 

9.8 

0.8 

1.7 

-1.9 

CRQ 

90071621 

7.0 

-0.8 

0.0 

4.0 

EMT 

90071621 
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63.0 

-0.0 

-0.0 

VCV 
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76.0 

54.0 

-0.3 

0.9 

lUiR 

90071008 

70.0 

65.0 

-0.0 

-0.0 

ONT 
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75.0 

59.0 

2.0 

3.5 

PMD 
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1.7 

9.8 
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70.0 

0.4 

65.0 

-1.0 

-0.0 

2.9 

-0.0 

-0.1 

WJF 
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71.0 

49.0 

4.9 

-0.9 

LGB 
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71.0 

-0.7 

64.0 

-2.0 

4.3 

-4.3 

-2.5 

2.5 

PMD 
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1.7 

9.8 
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VCV 
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75.0 

53.0 

-0.2 

1.0 
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BUR  90071809 

70.0 

65.0 

1 

o 

o 

O 

o 

1 

ONT  90071809 

75.0 

60.0 

5.9 

-1.0 

NZJ  90071809 
WJF  90071809 

70.0 

0.1 

62.0 

2.0 

-0.0 

4.0 

-0.0 

-0.8 

SDB  90071809 
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-5.4 
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-1.0 
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'  o 

3.0 

-3.0 
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BUR  90071810 

70.0 

65.0 

-0.0 

-0.0 

ONT  90071810 
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70.0 
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-0.0 

-0.0 
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72.0 
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-0.0 

-0.0 

NZJ  90071811 

69.0 
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63.0 

2.0 

-0.0 

6.5 
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WJF  90071811 

77.0 
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1.9 
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69.0 

0.3 
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BUR  90071812 
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ONT  90071812 
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-0.0 

NZJ  90071812 

69.0 
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63.0 

2.0 

-0.0 

6.5 
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3.7 

WJF  90071812 

77.0 

52.0 

10.8 

1.9 

SDB  90071812 

69.0 

50.0 

6.1 

-5.1 

LAX  90071813 

69.0 

0.3 

65.0 

-2.0 

5.2 

-5.2 

3.0 

-3.0 

RIV  90071813 

72.0 

63.0 

-0.0 

-0.0 

SBD  90071813 

71.0 

63.0 

-0.0 

-0.0 

VCV  90071813 

73.0 
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-3.1 

2.6 

BUR  90071813 

70.0 

65.0 

-0.0 

-0.0 

ONT  90071813 

72.0 

59.0 

-0.0 

-0.0 

NZJ  90071813 

69.0 

-0.3 

63.0 

2.0 

-0.0 

6.5 

-0.0 

3.7 

WJF  90071813 

77.0 

52.0 

10.8 

1.9 

SDB  90071813 

69.0 

50.0 

6.1 

-5.1 

Obs  Mean 

; 

76.6 

60.8 

4.5 

2.5 

Error  Mean 

: 

-0.3 

0.4 

0.4 

0.3 

Abs  Error 

Mean : 

2.0 

1.9 

2.9 

3.5 
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*  Station  Summary  * 

***A******************)lf*5k*************** 


St alien  :  SDB 


Obs  Mean 

72.8 

52.8 

4.1 

0.9 

Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Abs  Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Station  :  MWS 

Obs  Mean 

72.2 

51.2 

1.3 

1.9 

Error  Mean 

3.9 

9.6 

1.4 

1.3 

Abs  Error  Mean 

3.9 

9.6 

1.7 

2.1 

Station  :  WJF 

Obs  Mean 

82.6 

54.1 

9.7 

9.3 

Error  Mean 

-1.6 

3.0 

-0.4 

-2.3 

Abs  Error  Mean 

2.0 

3.1 

1.6 

4.2 

Station  :  SMO 

Obs  Mean 

72.6 

65.0 

6.3 

6.1 

Error  Mean 

1.5 

-1.6 

-0.9 

-3.4 

Abs  Error  Mean 

1.7 

1.7 

2.2 

3.7 

Station  :  ONT 

Obs  Mean 

79.8 

59.8 

6.4 

1.6 

Error  Mean 

-5.3 

2.4 

-2.0 

-0.2 

Abs  Error  Mean 

5.3 

2.4 

3.2 

2.0 

Station  :  EDW 

Obs  Mean 

86.5 

59.1 

7.0 

5.3 

Error  Mean 

-999.9 

-999.9 

1.6 

4.1 

Abs  Error  Mean 

-999.9 

-999.9 

2.7 

4.8 

Station  :  VNY 

Obs  Mean 

77.7 

62.8 

-1.7 

2.0 

Error  Mean 

-2.6 

1.1 

2.6 

2.4 

Abs  Error  Mean 

2.7 

1.2 

3.2 

2.9 

Station  :  TOA 

Obs  Mean 

72.2 

-999.9 

4.9 

-2.2 

Error  Mean 

2.2 

-999.9 

1.0 

4.8 

Abs  Error  Mean 

2.3 

-999.9 

4.0 

5.4 

Station  :  SNA 

Obs  Mean 

74.8 

64.8 

3.1 

6.6 

Error  Mean 

-0.2 

-2.0 

1.0 

-5.4 
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Abs  Error  Mean; 

1.4 

2.0 

2.5 

5.5 

Station  :  RAL  _  _ 

Obs  Mean  ; 

-999.9 

-999.9 

10.3 

-J.Z 

Error  Mean  : 

-999.9 

-999.9 

-3.5 

4 . 4 

A-bs  Error  Mean: 

-999.9 

-999.9 

3.7 

4 . 4 

Station  :  POC 

Obs  Mean  : 

-999.9 

-999.9 

6.6 

-0.2 

Error  Mean  : 

-999.9 

-999.9 

0.4 

3.4 

Abs  Error  Mean: 

-999.9 

-999.9 

1.9 

3 . 5 

Station  :  PMD 

Obs  Mean  : 

-999.9 

-999.9 

8.8 

8.9 

Error  Mean  : 

-999.9 

-999.9 

-0.9 

-0.4 

Abs  Error  Mean: 

-999.9 

-999.9 

2.3 

3.4 

Station  :  HHR 

Obs  Mean 

-999.9 

-999.9 

7.1 

1.0 

Error  Mean 

-999.9 

-999.9 

-0.2 

2.0 

Abs  Error  Mean 

-999.9 

-999.9 

1.8 

2.2 

Station  :  FUL 

Obs  Mean 

76.5 

62.7 

2.9 

4.8 

Error  Mean 

-1.5 

0.4 

2.2 

-2.7 

Abs  Error  Mean 

1.9 

1.1 

3.1 

3.0 

Station  :  EMT 

Obs  Mean 

-999.9 

-999.9 

2.4 

4.3 

Error  Mean 

-999.9 

-999.9 

0.6 

-3.4 

Abs  Error  Mean 

-999.9 

-999.9 

1.8 

3.9 

Station  :  CRQ 

Obs  Mean 

:  -999.9 

-999.9 

4.0 

0.6 

Error  Mean 

:  -999.9 

-999.9 

2.9 

1.2 

Abs  Error  Mean 

:  -999.9 

-999.9 

5.3 

2.4 

Station  :  CNO 

Obs  Mean 

:  -999.9 

-999.9 

7.6 

3.8 

Error  Mean 

:  -999.9 

-999.9 

-0.0 

-2.6 

Abs  Error  Mean 

:  -999.9 

-999.9 

2.6 

3.4 

Slat  ion  :  BUR 

Obs  Mean 

:  74.7 

64.2 

-1.4 

3.6 

Error  Mean 

:  -0.8 

-1.1 

3.3 

-2.6 

Abs  Error  Mear 

:  1.5 

1.4 

3.8 

2.7 

-58- 


station  :  VCV 


Obs  Mean 

80.8 

54.1 

CM 

O 

3.7 

Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Abs  Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Station  :  SBD 

Obs  Mean 

79.0 

63.8 

2.5 

0.1 

Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Abs  Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Station  :  RIV 

Obs  Mean 

78.7 

63.8 

4.3 

-1.8 

Error  Mean 

0.7 

-0.4 

1.5 

2.2 

Abs  Error  Mean 

1.2 

1.7 

2.2 

2.2 

Station  :  NZJ 

Obs  Mean 

72.8 

62.7 

2.8 

0.6 

Error  Mean 

1.0 

1.5 

2.1 

3.0 

Abs  Error  Mean 

1.6 

1.6 

3.5 

3.4 

Station  :  NUC 

Obs  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Abs  Error  Mean 

-999.9 

-999.9 

-999.9 

-999.9 

Station  :  LGB 

Obs  Mean 

74.0 

64.5 

4.6 

-0.5 

Error  Mean 

-0.7 

-1.4 

-0.7 

2.5 

Abs  Error  Mean 

1.5 

1.6 

2.8 

5.0 

Station  ;  LAX 

Obs  Mean 

71.0 

63.9 

7.4 

2.8 

Error  Mean 

1.0 

-0.1 

-3.6 

-1.6 

Abs  Error  Mean 

1.4 

1.1 

4.0 

2.0 
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